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1.  INTRODUCTION  .AND  SUMMARY, 


Let  (X  Y.) ,  i  =  l,  ...,  n  be  independent  identically  distributed  pairs  of  positive 
random  variables  with  distribution  functions  F(t)  *  P(X^  s  t)  and  G(t)  =  P(Y^  s  t) .  Suppose 
that  X^  and  Y^  are  independent  for  each  i.  In  the  random  censorship  model  the  X's 
represent  survival  times  and  the  Y's  represent  censoring  times;  the  observations  consist 
of  (Z^ ,  5^) ,  i«l,  n,  where  =  min^Y^  and  6.^  =  I(Xi  s  Yi)  . 


The  most  commonly  used  estimate  of  F  is  the  Kaplan-Meier  estimator  defined  by 

.  y*(i) 

■ 1  -  n  j  ci.D 


Z(i)st 


where  Z,n  S  ...  SZ,  denote  the  ordered  values  of  Z.  .  Z.,,  ....  Z  ,  and  <5 , . ,  is 

(1)  C  2j  (n)  1*  ’  n'  (l) 

the  5  corresponding  to  Z^  .  Note  that  this  version  of  the  Kaplan-Meier  estimator  is 
strictly  positive  on  [Z^,«)  if  the  last  observation  is  censored  and  is  therefore 
slightly  different  from  the  version  which  is  always  0  on  [Z^,®),  originally  proposed 
by  Kaplan  and  Meier  (1958);  see  Wellner  (1985)  for  a  discussion. 

An  important  problem  in  survival  analysis  is  to  estimate  the  pth  quantile  of  F 


defined  by  £p*inf(t:  F(t)  ip}.  The  natural  estimate  of  Cp  is 

£p*  inf (t:  F(t)  ip}  . 


(1.2) 


The  purpose  of  this  paper  is  three- fold.  First,  in  Section  2  we  show  how  the 
theory  of  counting  processes,  martingales  and  stochastic  integration  can  be  used  to 
obtain  in  a  simple  way  the  following  result. 


THEOREM  1.  Suppose  that  the  following  conditions  hold  for  p  e  (0,1). 

(i)  F  is  continuous  and  G(S  )  <  1. 

P 

(ii)  F  has  a  positive  derivative  at  5  | 


P-FU) 

Z  *  z  *  - +  R 

’P  P  F'Up)  n,p  ’ 


(1.3) 


where  n  '  R  -  )  in  probability  as  n  -■»  . 
a  i  P 


Dl.' 

w 


Since  n  (F(Cp)  -  p)  converges  in  distribution  to  a  normal  random  variable  with  mean  0 


and  variance 


u  -  p)  p- 

J„  (1 


P  _ d  F  (u) _ 

J0  (1-  G(u))2(l-  F(u)) 


(1.4) 


y  * 

(see,  e.g.  Gill,  1983)  heorem  1  immediately  yields  that  n^(5^  -  £p)  is  asymptotically 
normal  with  mean  0  and  variance 


1  -  P 

FTF) 


P  '  0 


dF(u) _ 

(l-G(u) )  2 (l-F(u)  ) 


(1.5) 


We  remark  that  (1.3)  is  neither  new  nor  the  strongest  version  of  this  result.  Aly, 

Csorgo  and  Horvath  (1985)  obtained  under  more  restrictive  conditions  than  those  given 

in  Theorem  1,  a  uniform  (in  p)  rate  of  convergence  for  R  of  0(n~  4(logn)  2  (log  log  n)  4), 

n  ,p 

as  n  -*■■*>  (see  Proposition  3.1  of  Section  3  in  the  present  paper  for  a  statement  of  this 

result).  Aly  et  al.  (1985)  assume  in  addition  to  the  conditions  given  in  Theorem  1 

that  G  is  continuous  and  that  F  is  twice  differentiable  at  £p  in  order  to  obtain  for 

fixed  p  their  rate  of  convergence  for  R  .  Therefore,  Theorem  1  is  of  interest 

n,p 

because  of  its  simple  proof  and  because  its  assumptions  are  less  restrictive  than  those 
of  Aly  et  al .  (1985) . 

Our  second  purpose  is  to  prove  the  asymptotic  independence  of  the  estimates  of 

quantiles  in  the  competing  risks  problem;  this  is  done  in  Section  3.  The  competing 

risks  problem  arises  in  medical  follow-up  studies  and  in  industrial  life  testing  when 

the  failure  of  an  individual  can  be  classified  into  one  of  m  mutually  exclusive  classes, 

referred  to  as  "risks";  see  Gail  (1975).  Consider  a  sample  of  n  individuals  each  of 

whom  can  incur  a  failure  by  any  one  of  ra  risks.  Let  ,  i=l,  n,  j  =  l,  m 

denote  the  lifelength  of  individual  i  if  risk  j  is  the  sole  cause  of  failure.  Let  F^ 

denote  the  distribution  of  X...  Assume  that  X..,  ....  X.  are  independent,  and  that 

ij  ii  lm  r 

the  vectors  (X^ , . . .  ,X^m)  i=l,  ...,  n  are  indeper^en"  .  Since  the  life  of  individual 
i  is  terminated  by  the  earliest  occurring  cause  of  failure,  one  observes  only  the 
sequence  j  =  1 , . . .  ,m,  i=l , . . .  ,n} ,  where  Zi  =  min(Xi;L , . . .  ,Xim)  and  6i  j  =  I  (zi  =  xi  j )  - 


Thus,  <5^  indicates  whether  risk  j  was  the  cause  of  death  (6^  =  1)  or  not  (o^j  =  0)  . 


Note  that  for  each  j,  we  can  estimate  by 

F.Ct)  .  1  -  n  f-i  f(i,j  , 

3  z(i)st '  11-1+1  • 

where  Z^  s  s  •••  5  are  the  ordered  values  of  Z^,  Z^, 
the  corresponding  values  of  <5^,  &2y  *  $nj* 

The  pth  quantile  of  F^  may  be  estimated  by 

=  inf(t :  F.  (t)>p)  . 

J  >r  J 


(1.6) 


Z  and  <5 , .  .  are 
n  UJl 


(1.7) 


We  combine  the  Bahadur  representation  of  quantiles  given  by  Aly,  et  al.  (1985)  and  a 
weak  convergence  result  of  Aalen  (1976)  to  show  in  a  simple  argument  that  when  suitably 

A  A 

normalized,  the  vector  (51  ,...,£  )  converges  in  distribution  to  a  vector  of  inde- 

l ,  p  m ,  p 

pendent  Gaussian  processes. 

Finally,  in  Section  4  we  use  the  ideas  in  the  proof  of  Theorem  1  to  study  the 
quantile  estimates  that  arise  in  the  reliability  model  introduced  by  Doss,  Freitag, 
and  Proschan  (1986). 

2.  A  REPRESENTATION  FOR  QUANTILES  OF  THE  KAPLAN-MEIER  ESTIMATOR 

An  account  of  the  theory  of  counting  processes,  martingales,  and  stochastic 
integration  used  in  this  section  is  given  in  Chapter  18  of  Lipster  and  Shiryayev  (1978) 
and  in  Chapter  2  of  Gill  (1980) .  A  very  accessible  exposition  of  this  theory  is  pro¬ 
vided  in  Anderson  and  Borgan  (1984).  Throughout  the  paper  we  adopt  the  convention 
that  ^  =0;  xAy  denotes  min(x,y),  and  if  H  is  a  distribution  function  then  H(t)  =  1  -  H(t) 

The  following  processes  are  needed  to  give  (F(t)  -  F(t  A  Z,  ,))/F(tAZ,  ,)  a  stochast 

(.nj  (nj 

integral  representation  which  we  will  use  to  prove  our  results. 

J(t)  =  I(Zr  .  >t)  (2.1) 


V(t)  =  l  I  (Z  >  t) 
i=l 


(2.2) 


1 


N(t)  =  l  HZ  <  t,  6  =  1) 
i  =  l  1 


(2.3) 


&  . .  . 

/  /  /  /  / 


.  :• .  - 


i 


\ 

♦ 

I 


A(t)  =  / 


0  F(s) 


dF(s) 


(2.4) 


M(t)  =  N(t)  -  A(t) 


(2.5) 


Our  proof  of  Theorem  1  will  require  the  following  two  results  taken  from  the 


literature. 


PROPOSITION  2.1  Suppose  F  is  continuous.  Then 


—  :  F(tAZw} .  f  hizUM  dM(s) 

%,t]  p(s> v(s) 


(2.6) 


for  tiO,  where  P(s-)  =  lim  1  -  F(s) ,  and 

u+s 


F(s-)  J(s) 


(0,t] 


PCs)  V(s) 


dM(s) 


2 

n*  PCs-) ,J(s)  d  A(s) ,  F  ;  t  €  [0, 

F(s)  V (s)  t 


(0,t] 


is  a  martingale,  where  F  is  the  completion  of  the  o-field  generated  by 


a(I(Z^  s  s,  <5^  =  1) ,  I(Z^  s  s,  6^  =  0) ;  lsisn,  s  s  t)  and  T  >  0  is  such  that  F(T)  <  1. 


PROOF.  The  identity  (2.6)  is  given  in  Lemma  3.2.1  Part  (iv)  of  Gill  (1980),  and  also 


in  Theorem  3.1  of  Aalen  and  Johansen  (1978). 


To  prove  Part  (ii)  ,  we  note  that  it  follows  from  Corollary  3.1.1  of  Gill  (1980) 


that  { (M(t) ,  Ft )  ;  te  [0,T]}  is  a  martingale  with  quadratic  variation  process  <M>(t)  =  A(t) 
Let  f(s)  =  (F(s-)  J(s))/(F(s)  V(s)).  Since  f  is  predictable ,  the  integral 


/  f(s)  dM(s)  is  a  martingale  with  quadratic  variation  process  J  f2(s)  d<M>(s), 

(o , t]  (0,t] 


and  Part  (ii)  follows  (see  p.  268  of  Lipster  and  Shiryayev,  1978  or  p.  10  of  Gill,  1980; 


cf.  also  the  paragraph  preceding  (**.1.11)  page  55  of  Gill,  1980). 


PROPOSITION  2.2  Let  (V  }  and  (W  }  be  two  sequences  of  random  variables  satisfying 


n  n=l 


the  following  conditions. 


k 

$ 

i 


(i)  For  all  a>0  there  is  a  X  (depending  on  o)  such  that  P(  |Wn  |  >A)  <  a  (2.7) 
for  all  n. 


(ii)  For  all  k  and  e  >  0, 


lim  P(^n  Sk,  W  2  k  +  e)  =  0  and  lim  P(Vn  2  k  +  e,  £  k)  =  0 
n-+®  n-«° 


(2.8) 


Then  vn  “  wn  "*■  0  in  probability  as  n-*-*. 


PROOF:  See  Ghosh  (1971). 


PROOF  OF  THEOREM  1.  We  note  that  the  assumption  F'(C^)  >  0  implies  that  is  the  unique 
number  satisfying  F(£  )  =p.  For  any  t  i  0  and  any  n  define 


Z,  n  3  n%CF^  -  F(E  ♦  3/F'Uj  ,  and 

t,n  P  n'*  P  n't  P 


Note  that 


where 


Wn«n^(p-  F(5p))/F-(S  )  . 


<n‘,Clp.«p)st>.<psFC5p.-Sf)>.  Ut>nStn}  , 


(2.9) 


(2.10) 


by  definition  of  the  derivative.  In  (2.9)  the  first  equality  follows  from  the  definition 


of  given  by  (1.2).  Suppose  we  can  show  for  fixed  t20  that 


Z..  -  W  -*■  0  in  probability. 

t,n  n  r 


(2.11) 


Letting  Vn  =  n^(£p  -  5p) »  it  would  then  follow  from  (2.9)  and  (2.10)  that  (2.8)  holds  for 
the  sequences  {V^}  and  {W^} .  Since  W^  is  asymptotically  normal  (see  e.g.  Theorem  4.2.2 
of  Gill,  1980)  condition  (2.7)  is  satisfied,  and  the  theorem  follows  from  Proposition 


It  remains  to  prove  (2.11).  Algebraic  manipulations  yield  that 


Z  -W  =  -  F(5  +-V)-F(C) 

t,n  n  F  (C  )  [  P  ^  PJJ 


F(CP) 


(2.12) 


*  m  ♦ 4- )  FUy  ^ 3 .  Fiv 


F^v 


P 


F(C_  ♦  -V)  F(£  ) 

p  ^  p 


The  first  product  term  on  the  right  side  of  (2.12)  converges  to  0  in  probability  by  the 

consistency  of  the  Kaplan-Meier  estimate  at  5p  (see  e.g.  Foldes,  RejtcJ,  and  Winter, 

1980)  and  the  definition  of  the  derivative. 

We  now  proceed  to  show  that  the  second  product  term  on  the  right  side  of  (2.12) 

converges  to  0  in  probability.  For  fixed  t,  there  exists  N  (depending  on  t)  such 

that  for  n>N,  F(C  +-£-)<  1.  Thus,  for  n^N,  by  Proposition  2.1,  the  second  product 
P 

term  on  the  right  side  of  (2.12)  is  equal  to 

J'i  +  T  Tfc'l 


F'(Cp)  F(CP  + 


t  ,  F(s-)  J< 

2*  .£  ♦  4-1  P(S)  V( 


dM(s) 


(2.13) 


- 

F(«H 

-V) 

P(5  + 

P 

P 

P(£ 

►  4-  ) 

P((*L 

l  P 

n^ 

P 

'F(y  Z(n)) 


It  is  clear  that  under  the  assumptions  of  the  theorem,  the  probability  that  the  quantity 

inside  the  brackets  in  (2.13)  is  not  0  converges  to  0  (exponentially  fast)  as  n 

We  now  show  that  the  first  term  in  (2.13)  converges  to  0  in  probability.  Let 

e  >  0  be  such  that  F(£  +  e)  <  1,  and  G(?  +  e)  <  1.  By  the  strong  uniform  consistency  of 

P  P 

the  Kaplan-Meier  estimate  (see  Foldes,  Rejto,  and  Winter,  1980), 


sup  |  -1-0 

Os  ss  (  +  e  F(s)  j 
P 

It  follows  from  the  Glivenko-Cantel li  Theorem  and  (2.14)  that 

I —  “T  2 


um  sup  n"2  dA(s) 

**•  F(s)  J 


P  P 


d  F(s) 


H(s)  F ( s J 


(2.14) 


(2.15) 


where  H(s)  i  P(Ij>s).  Note  that  H(s)  is  bounded  away  from  0  on  [%  ,  5  +e). 


mm® 


A  V,/ 

%  W*-*  m  w  J  *•  - 


■vvN)/!vj 


Since  e  was  arbitrary,  (2.15)  implies  that 


.vVL“ 


%  F(s-)  J(s) 
F(s)  V(s) 


dA(s)  -*■  0  in  probability.  (2.161 


We  now  apply  Lenglart's  inequality  (Lenglart,  1977;  for  a  statement,  see  Theorem  2.4.2 
of  Gill,  1980)  to  the  martingale  of  Proposition  2.1,  Part  (ii)  to  conclude  that 


?(S~j  d  M(s) 

Vp*T]  F(s) 

n 


in  probability. 


This  completes  the  proof  of  the  theorem.  ■ 

The  proof  of  Theorem  1  generalizes  the  method  used  by  Ghosh  (1971)  to  the  case 
of  censored  data. 

v  - 

To  estimate  the  expression  for  the  asymptotic  variance  of  nJ(£p-  given  by  (1.5), 

we  may  estimate  the  term  (1  -  p) 2/ (F '(£  ))2  and  the  integral  term  separately.  Under 

P 

the  conditions  of  Theorem  1,  the  proposition  in  Section  2  of  Hall  and  Wellner  (1980) 

r  5(i) 

implies  that  n  )  — : — rr7 — rr  is  a  consistent  estimate  of  the  integral  in 

i:Z,.Se  (n-i+1)  (n-i) 


(1.5). 


x,ZCi)  5p 

If  we  assume  in  addition  to  the  conditions  of  Theorem  1  that  F  is  continuously 


differentiable  in  a  neighborhood  of  5^,  then  Theorem  4.1.2  of  Ramlau-Hansen  (1983) 

implies  that  F'(C  )/F(£  )  may  be  consistently  estimated  by  a  suitable  kernel  type  of 

estimator;  one  such  estimator  is  described  in  Section  5  of  Ramlau-Hansen  (1983). 

Note  that  since  E,  is  unknown,  it  must  be  replaced  bv  Z  in  the  estimate  of 

p  P 

*»  J-'  a 

F'(C  )/F(C  ).  Thus,  the  asymptotic  variance  of  nJ(£  -Z  )  can  be  consistently  esti- 
P  P  P  P 

mated,  enabling  the  construction  of  asymptotic  confidence  intervals  for  £  . 

We  remark  that  for  small  samples  the  confidence  intervals  for  Z  based  on  the 

above  method  are  not  aopropriate.  Several  methods  have  been  proposed  for  setting  srr^ll 

sample  confidence  limits  for  C,^;  see  Slud,  Byar,  and  Green  ■;  1984  '  ,  and  the  references 

2 


cited  therein. 


3. 


WEAK  CONVERGENCE  OF  THE  QUANTILES  OF  THE  KAPLAN-MEIER  ESTIMATORS  IN  THE  COMPETING 
RISKS  PROBLEM. 

In  this  section  we  use  the  representation  i . (p)  =  £ . (p) + (p-F . (5. (p) )) /F?(Z- (p) ) +  R-  _Cp) 

J  J  J  J  J  J  J  >1L 

to  obtain  in  a  simple  way  the  asymptotic  distribution  of  (£^(p),  . Sm(p))*  In  the 

sequel  we  use  the  notation  £.(p)  rather  than  t-  _  when  the  quantiles  of  F.  are  viewed  as 

J  1  >P  J 

a  process;  similarly  for  £.(p).  Let  G.  denote  the  distribution  of  min  X,  .  The  following 

^  J  k*j 

proposition  due  to  Aly  et  al.  (1985)  gives  the  strongest  known  result  concerning  the 

remainder  process  R.  (p) . 

F  J,n  K 

PROPOSITION  5.1.  Suppose  that  the  following  conditions  hold. 

(i)  F.  and  G.  are  continuous. 

J  J 

(ii)  F^  is  twice  differentiable  on  (t^.T^),  where  t^  =  sup{x:  F  (x)  =0 } 
and  Tj  =  inf{x:  F^.(x)  =  l}. 

(iii)  fj(t)=Fjlt)  is  positive  on  (tj,Tj). 

(iv)  For  some  P*  e  [0,1]  we  have 

sup  p|fru.Cp))|/(f-U.(p))2  <- 
OSpSp*  J  J  J  J 

(v)  0  <  lim  f .  (t)  <  °°  . 

t+t.  J 
J 

Then,  if  pQ  <  p*  and  G^  (£ ^  (p0) )  <  1 ,  we  have  as  n+" 

sup  |R.  (p)  |  =  0(n'/4(log  n f  (loglogn)4)  a.s. 

0*P*P0  J'n 

PROOF.  See  Corollary  4.1  of  Aly  et  al.  (1985). 

For  s  >  0,  let  D[0,s]  be  the  space  of  all  real  valued  functions  defined  on  [0,sj 
that  are  right  continuous  and  have  left  limits,  with  the  Skorohod  metric  topology,  and 
let  DmC0,s]  denote  the  product  metric  space. 

THEOREM  2.  Suppose  that  the  following  conditions  hold  for  each  j  =  1,  ...,  m. 

(i)  Fj  is  continuous. 

(ii)  Fj  is  twice  continuously  differentiable  and  has  a  positive  density  on 

(t. ,T.) . 

J  J 


’&&& . . 


p|fj(5j(p))l 

(iii)  For  some  p*  £  [0,1]  we  have  sup  - — -  <  00 

Ospsp*  (f.(€  (p))T 

(iv)  0  <  lim  f . (t)  <  <*>  . 


(v)  P0  £  [0,p*)  is  such  that  max  CpQ^ )  <  1- 


Then,  as  n  -*■  «, 


weakly  in  Dm[0,pg],  where  U^,  ...,  are  independent  mean  0  Gaussian  processes  with 
covariance  structure  given  by 


Covflh  (p1) ,  Uj  Cp2) )  « 


(l-Pl)(l-p2)  (  dw _ 

WP133  fjCS(P23  3  .  TTFk(C.(w))  (1-w) 

0  k=l  J 


for  0  s  Pl  s  p2  s  pQ 


PROOF.  It  follows  from  Theorem  10.1  of  Aalen  [1976)  that 


P-  F1Cg1CP)) 
f1cc1(p)) 


P-FmCCm(P)3 

fm^m(p)) 


[3.1) 


(U. [p) ,  ...,  U  (p) )  weakly  in  D  CO ,pn) 


The  theorem  follows  from  (3.1)  and  Proposition  3.1. 

Theorem  2  is  of  statistical  interest  because  it  enables  the  construction  of 


simultaneous  confidence  intervals  for  C^(p^),  ...,  5mCPm)* 

We  remark  that  for  fixed  p£  (0,1),  if  F^  is  continuous,  Fj*  exists  and  is  posi¬ 
tive  at  C.  ,  and  max  F.  (£  .  )  <  1  for  1  =  1,  . . .  ,  m,  then  as  n+® 

3  »P  k  k  J.P 

u  * 

n  (Cx  -  q  ,  •  •  •  »  5m  p  -  p)  -*•  (u^p)  ,  . . .  ,  Um(p))  in  distribution.  This  is  easily 
seen  by  substituting  Theorem  1  for  Proposition  3.1  in  the  proof  of  Theorem  2. 


9 


& 


S 


4.  ESTIMATING  THE  QUANTILES  OF  A  COHERENT  SYSTEM. 


Consider  a  coherent  structure  of  m  independent  components.  For  the  sake  of 
concreteness,  we  keep  in  mind  the  specific  example  given  by  Figure  1  below. 

I 


Figure  1,  A  coherent  system  of  3  components. 
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Let  F  denote  the  distribution  function  of  the  lifelength  of  the  structure.  Doss, 
Freitag,  and  Proschan  (1986)  (subsequently  referred  to  as  DFP)  study  the  problem  of 
estimating  F  under  the  following  model.  A  sample  of  n  systems,  each  with  the  same  struc¬ 
ture  is  available  for  testing.  Each  system  is  continuously  observed  until  it  fails. 

For  every  component  in  each  system,  either  a  failure  time  or  a  censoring  time  is  recorded. 
A  failure  time  is  recorded  if  the  component  fails  before  or  at  the  time  of  system  failure. 
A  censoring  time  is  recorded  if  the  component  is  still  functioning  at  the  time  of  system 
failure. 

DFP  propose  the  following  procedure.  Let  F^ ,  ...,  Fm  be  the  distribution  functions 
of  the  lifelengths  of  the  m  components.  Write 


) 


P(t)  •  h(F1(t),  ....  Fm(t))  fort>0,  (4.1) 

where  h:  [0,l]m  -*■  [0,1]  is  the  reliability  function  (see  Chapter  2  of  Barlow  and 
Proschan,  1981,  for  details  concerning  reliability  functions).  In  the  example  given  by 
Figure  1,  it  is  easy  to  check  that 

F(t)  =  F1(t)  [  1  -  (1-F2(t))(l-F,(t)] 


so  that 


h(u1,u2,u3)  =  UjCl-  (1-^)  (1-u,)  ]  for  ux,  ,  u3  e  [0,1] 


Let  F 


...,  F^  be  the  Kaplan-Meier  estimates  of  F ^ ,  ...,  F  .  To  estimate  F,  DFP  propose 


F(t)  =  h(F1(t) ,...,Fm(t)) 


(4.2) 


10 


US  mV  m  ^  jt  » 


V/-" 


as  an  alternative  to  the  naive  estimate  given  by  the  proportion  of  systems  still 
functioning  at  time  t. 


The  purpose  of  this  section  is  to  study  the  estimates  of  quantiles 

t( p)  =  F ~ 1  Cp)  and  l  =  F-1  (p)  (0<p<l)  .  (4.1) 

J  J  P 

Before  doing  so  we  need  to  discuss  the  estimates  and  F. 

The  weak  convergence  of  the  Kaplan-Meier  estimator  to  a  Gaussian  process  has  been 
well-established  in  the  literature  (Breslow  and  Crowley,  1974;  Aalen,  1976;  Gill,  1983) 
under  the  assumption  that  the  lifelengths  and  the  censoring  variables  are  independent. 

In  our  situation  the  component  lifelengths  are  censored  by  the  system  life  length,  and 
the  independence  condition  is  clearly  violated.  We  can,  however,  redefine  the  censoring 


variables  to  bypass  this  difficulty.  This  is  easiest  to  explain  in  terms  of  the  example 
given  by  Figure  1.  Let  X^  =  lifelength  of  component  j.  Consider  Component  1.  Clearly, 
is  censored  by  =  X2  v  X^,  which  is  independent  of  X^.  Similarly,  X7  is  censored 
by  Y7 =  X^ ,  and  X,  by  Y^  3  X1 .  The  construction  is  general:  for  an  arbitrary  system,  X^ 
is  censored  by  Y ^  =  lifelength  of  system  if  X^.  is  replaced  by  DFP  show  that 
(i)  Xj  and  Y^  are  independent. 


(ii)  The  distribution  function  of  Yj  is 


G.(t)  =  1  -  h(F. (t) 


,...,  F.^Ct),  1,  Fj  +  1(t),  ...,  Fm(t)) 


(4.4) 


The  main  result  of  DFP  can  now  be  stated. 


PROPOSITION  4.1.  Suppose  that  F^,  — ,  F  are  continuous,  and  let  T  satisfy  max  F.(T)  <  1, 

1  UjSm  1 


Then  as  n  -*■  «  , 

n*^l"Fl’  F2-F2,  ...,  Fm-Fm)  -  (Wx,  W2 ,  . ...  WJ 

weakly  in  D  [0,T],  where  W  . ..,  W^  are  independent  mean  0  Gaussian  processes.  The 
covariance  structure  of  is  given  by 

C1  dF.(u) 

(i)  Cov(W  (t  )  ,W  .(t9))  =  F  .(t  )F.  (t  )  /  — -  for  0<t  <t,<T. 

J  1  1  2  J  1  J  2  0  GjluHFj(u))2  1  - 


n  J(f  -  F)  -*■  W  weakly  in  D[0,T], 


where  W  is  a  mean  0  Gaussian  process  with  covariance  structure  given  by 

m 


Cov(W(t 

where 


t,  dF.(u) 

1),W(t2))  =  ^h.(t1)h.(t2)  F.Ct^F.Ct^  /  7—7 L=r~  .2 


G.(u)(F  (u)) 


3h 


hj  ^  =  3u.  ^ul»*  * • »unp 
J 


(4.5) 


(u^,...,um)  a(Fj(t),...,Fm(t))  • 

The  model  being  considered  is  a  generalization  of  the  competing  risks  model,  which 
corresponds  to  a  series  system  of  m  components.  The  following  result  generalizes 
Theorem  2  to  this  model. 


PROPOSITION  4.2.  Let  5^.  (p)  be  defined  by  (4.3),  and  suppose  that  the  conditions  of 


Theorem  2  are  satisfied.  Then  as  n  +  ® 


n2(^l  '  5i»  ^Ul»  *•*»  UnP 


weakly  in  Dm[0,p0],  where  U^,  ....  are  independent  mean  zero  Gaussian  processes. 


The  covariance  structure  of  IK  is  given  by 


Cov(lL  (px) ,  U^(p2))  = 


(1-PX) Cl-P23 


fjCSjCPi))  fj(5j(p2)3 


dw 


Gj (Cj (w) )  (1-w)2 


for  OSPlsp2Sp0 


where  G^  is  given  by  (4.4). 


PROOF.  Substitute  Part  (i)  of  Proposition  4.1  for  Theorem  10.1  of  Aalen  (1976)  in  the 
proof  of  Theorem  2.  — 


The  main  result  of  this  section  gives  a  Bahadur  representation  for  analogous 


to  (1.5),  and  is  useful  in  obtaining  the  asymptotic  distribution  of  C  . 


THEOREM  3. 


Let  0<p<  1  and  assume  the  following  conditions. 


Then, 


(i)  F^,  ....  Fffl  are  continuous,  with  max  F.(C  )  <1 

lsjsm  J  p 


(ii)  For  each  j,  F.  has  a  positive  derivative  at  £  . 

3  P 


P*F«J 

S  *  S  *  F'Up)  *  Vp  ’ 


where  n2  Rn  ■+  0  in  probability  as  n-*-®  . 


PROOF.  The  proof  we  give  is  a  straightforward  modification  of  the  proof  of  Theorem  1, 
Since  h  is  strictly  increasing  in  each  of  its  arguments  (see  Theorem  1.2  of  Barlow 
and  Proschan,  1981)  it  follows  that  is  the  unique  number  satisfying  F(C  )  = p. 

Furthermore,  h  is  differentiable  over  [0,l]m  (see  Lemma  2.1  of  DFP) ,  so  by  the  chain 
rule  F'(S  )  exists  and  is  given  by 


F'(V  •  1VV  Y(V 


where  h^  is  given  by  (4.5). 

For  any  t  2  0  and  any  n  define 


fv£)  •  /f-(5p)  “d 


--'cp 


We  have  (n^-s^st)  =  .here  t„  =  ^  [f«p  ♦  ^  )  .  p]/F.(  ^  - 1 

(cf.  equations  (2.9)  and  (2.10)).  As  in  the  proof  of  Theorem  1,  it  suffices  to  show 
that  Z'  n  -  W'  ■*  0  in  probability. 

Using  (4.2)  we  write 


Z'  -  W*  =  — - - 

t,n  a  F-rcp) 


h(_F(y^))  -h(FU  )) 


b(F(C  +  — p  )  -  h(F(5 _))  ,  (4.6: 


where  for  economy  of  notation  F(s)  and  F(s)  denote  the  vectors  (F^(s),  Fm C s j ) 

and  (F^(s),  ....  F^C s) ) ,  respectively.  We  now  apply  the  Mean  Value  Theorem  (see  for 

example  Apostol,  1964,  Theorem  6.17)  to  each  of  the  differences  inside  the  brackets  in 

(4.6):  there  exist  points  a  and  a  lying  on  the  line  segments  joining  F(?  +— r)  and 

"  ‘  P  n"* 

-  _  t 

F(£p) ,  and  F(£  ♦— )  and  F(£  ),  respectively,  such  that  the  right  side  of  (4.6)  is 

-  f  P  n't  P 

equal  to 


F  («p) 


Vh(a) 


'  H'l  -  7hta)-  favA’  -  *(v" 

J  L  n 


)  -  F (£  )  ,  (4.~) 


where  7h  denotes  the  gradient  of  h.  It  is  convenient  to  rewrite  (4.7)  as 


p'cy 


Vh(a) 


F(£  +  ~T"  )  -  F(£  ) 
-  P  rfi  J  -  P 


iVi 

H  n  • 


)  -  F(E  ) 


F-(C  ) 
P 


(Vh(a)  -  Vh(a)) 


> 

&V?5 --f(v] ;  • 


(4.3) 


Since  we  clearly  have  t  -*■  F(£^)  a.s.  and  a  -*■  F(£p),  the  continuity  of  Vh  (see  Lemma 
2.1  of  DFP)  and  the  differentiability  of  each  Fj  at  imply  that  the  second  term 
inside  the  braces  in  (4.3)  converges  to  0  in  probability.  The  first  term  inside  the 
braces  in  (4.8)  converges  to  0  in  probability  since  by  (2.11),  for  each  j 

'  > 

n'2  F.^*-T  )  -  F.  (C  )  -  F.(£  +4r)  -  F.  (  £  +-Tr)l  ^  0  in  probability, 

L  J  P  ni  J  p  -1  - ]  p  n  *  J  p  n  * 

and  the  components  of  7h  are  bounded  by  1  f see  Lemma  2.1  of  HFP'' .  This  completes 


the  proof  of  Theorem  3. 


Theorem  3  immediately  implies  the  asymptotic  normality  of 


COROLLARY  4.1.  Under  the  conditions  of  Theorem  3,  n^CC^  -  ■+■  NCO.a^p)),  where 


‘Cp)  ■  Fh?  ji  [hj‘V  WJ 


-  2  p  dF .  (u) 


G- (u) (F . (u) ) 


J 


To  construct  asymptotic  confidence  intervals  for  Cp  we  need  to  be  able  to  consistently 
estimate  the  quantity  a2(p).  An  extension  of  the  argument  given  at  the  end  of  Section  2 
shows  that  this  can  be  done  under  the  additional  assumption  that  for  each  j,  F.  is 
continuously  differentiable  in  a  neighborhood  of  5  .  Details  are  omitted. 
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